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I. INTRODUCTION 


Holographic duality (gauge/gravity correspondece) [iHSl is a particular kind of 
strong/weak duality, which gives us a powerful method to study strongly correlated sys¬ 
tems in condensed matter physics via the study of the black hole physics with one extra 
dimension. Especially, it provides a novel way to explore the strongly correlated fermionic 
system. 

The early works jd-S] on this subject have made efforts on understanding the mysteri¬ 
ous behaviors of the existing non-Fermi liquid which is named as “holographic non-Fermi 
liquid”. In detail, by studying the evolution of the bulk Dirac equation in the RN-AdS 
black hole with the ingoing boundary condition at the horizon, the fermionic correlator can 
be extracted at the AdS boundary holographically. This proposal attracts many interests 
in the related topics. Some endeavors have been made to generalize the investigations of 
the probe fermion to the general charged Gauss-Bonnet black hole, charged dilatonic black 
hole, charged Lifshitz black hole and so onjo-lS]. One also implemented the holographic 
non-relativistic fermionic hxed points by imposing the Lorentz violating boundary condi- 


tion, which results in an inhnite flat band in the dual boundary held theory 16l-ll8| . By 


introducing the interaction between Dirac held and gauge held, it was disclosed that the 
dual liquid can transform from the Fermi liquid to non-Fermi liquid and also to the Mott 
insulating phase 


19 


24j |. The properties of the fermionic spectra’ 


ehects, which is introduced in 25l-l27l|. has been explored in 28 


function including lattice 


29|. It shows some crucial 


characteristics, including periodic structure and Brillouin zones, which can be compared with 
the results of Angle-resolved photoemission spectroscopy (ARPES) or Scanning tunneling 
microscopy (STM), in real condense physical system. 

In real condensed matter system, another crucial role is impurity. The physical properties 
of materials will drastically change when impurity is added. A famous example is the 
cuprates, which exhibits Mott insulating phase, non-Fermi liquid phase and Fermi liquid 
phase as the doping levels change. Another well-known example is the Kondo effect that 
when the magnetic impurity is introduced, the resistance increases logarithmically with 
the decrease of the temperature. In holography, impurity can be introduced in Einstein- 
Maxwell-massless scalar gravity j^, M|, ^ which it is implemented by spatially dependent 
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sources for scalar field operators with being linear in coordinates^. By introducing an extra 


charged complex scalar field, the authors of 


M 


35| built a holographic superconductor with 


momentum relaxation. They found a new type of superconductor induced by the effective 
impurity a even at vanishing chemical potential p = 0. More interestingly, the effective 
impurity a in holography can induce a coherent/incoherent metal transition ^,[3^. The 
coherent metal has a standard Drude peak, which is usually described by Fermi liquid 
theory, while the incoherent metal follows non-Drude behavior, which complies with non- 
Fermi liquid theory. This motivates us to investigate the fermionic response with impurity 
in holography. 

Subsequently, we will shortly review the simple holographic model of momentum relax¬ 


ation 


30 


31| and discuss its near horizon geometry at any temperature in section II. Then 


the equation of motion and Green function on this background are derived. Afterwards, we 
exhibit our numerical results and discuss Fermi surface structure and phase diagram induced 
by the effective impurity parameter a. Finally, our results are summarized in section V. 


II. EINSTEIN-MAXWELL GRAVITY WITH MASSLESS SCALAR FIELD 
SOURCES 

We start with a brief review of the bulk Einstein-Maxwell gravity accompanied with 
massless scalar field sources and analyze its near horizon geometry. 


A. Background geometry with scalar field sources 

We combine the antion of free massless scalars together with the Einstein-Maxwell action 
in d -|- 1 dimension 


31|, 




d-l 


R-2A- - i 


A f d-‘x^2K, (1) 

JdM 


where A = —d{d — 1 )/( 2 L^) and the held strength = d^A^, — for a U{1) gauge 
held. The second term is the Gibbons-Hawking term where 7 is the induced metric on the 
boundary and K is the trace of the extrinsic curvature. 


^ Impurity is also implemented by disorder in holography. For instance, one can refer to 3^ [33 1. 
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For simplicity, we choose the gravitational coupling = IQiiGd+i and the AdS radius 
L as unity. The equations of motion can be easily obtained as 


d-l 


R 






R-2A- 


1 


d-l 




7=1 


7=1 


= V^^J = 0. (3) 


The equations of motions admit the following isotropic solution of various helds 


dr"^ 


d-l 


ds^ = —f(r)dt^ + —— + dx^dx\ 

f{r) ^ 


i=l 


r 2 ^ ^0 . 

^ ^ 2{d-2) r'^-2 




„d-2 


where 




mo = r“ 1 + 


dt, 


^2{d-2) ' 
T 7 = a/jX* = a6iix\ 


d — 2 jR 


a 


2{d — 1) Tq 2{d — 2 ) Tq 


Q^ = 


{d - 2)p2r2('^-2) 


(4) 

(5) 

( 6 ) 


2 (d-l) 

and o? = ■ *5* with the vector notation (cfi)/ = an and (5* • a^ = 'Ylii Oiuai^. We 

can see that due to the presence of the spatial dependence in the massless scalar held, the 
full solution is not isotropic and homogeneous, though the metric is isotropic. Note that the 


anisotropic solutions with only one scalar held has been studied in 40|, 1^. The parameter a 


mimics the strength of the translational symmetry, which means the momentum relaxation 
in the dual theory. Note that when a = 0, the solution goes back to the RN-AdS black hole 
dual to the held theory with translational invariance. 

The temperature and entropy density of the black hole are given by 


y _ firp) ^ ^ ^ 

dvr dvr 


(d — 2)2Q^ 


2drl 


drl'^ ^ 


s = dvrr, 


d-l 


(7) 

( 8 ) 


2 Note that it was addressed in [SlJ that this solution is highly similar to the solution proposed in massive 
gravity 37 Q- 

^ It is worthy pointing out that in order to consider the CFT dual to the anisotropic system, we have to 
focus on the dimension d ^ 3. Specially, in the background with dimension d = 1, a logarithmic divergence 
will appear and a cannot be introduced into the action. 
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The solution is fully characterized by T, /i and a. For a given chemical potential and 
temperature, we can determine the value of a. This model can be seen as the case of the 
long wavelength limit of Q-lattice model [2^. But an periodic lattice structure is absent 
because the scalar source is linearly dependent on the spatial coordinates. Therefore, it is 
more reasonable to treat a as impurity, as illustrated in j^, 361. 


B. Near horizon geometry 


In order to study the near horizon geometry at any temperature, we will introduce a 
length scale r* by parameterizing Q with dimension d — 1 as 


Then the temperature can be written as 

4:71 




d 


„d-l 


2drl \ d-2 * 


a 


2drl 


l-(- 


* \2d-2 


(9) 


( 10 ) 


We can see that r* is equal to the horizon radius tq at zero temperature and smaller than 
the radius ro at finite temperature. To work out the near horizon geometry, we consider the 
scaling limits 

r2 r2 

^ ^ ( 11 ) 


r — r* = cu—, ,ro —r*=a;- , 

C Co 


t = (jJ 


with a; —0, C and r hnite. L 2 is the curvature radius of AdS 2 geometry whose expression 
will be shown later. 

Under the above transformation, the leading term of the redshift factor / in (j4]) is 




( 12 ) 


with the curvature radius a^+{d- 2 )'^% ^ • Thus, metric near horizon becomes 


r 2 

ds^ = — 
^2 


e 

c, 


1 -TTT 1 dv"^ + 


1 - 

‘'0 


d^ 


d-l 


dx^dx^. 


(13) 


2=1 


and the temperature is T = 2 ^- The zero temperature limit will be reached as Co 00 so 
that the formulation of metric can be deduced from Eq. (lT3ll as 

e 


ds^ = ^ {—dr'^ + dC) + rl ^ dx^dx\ 


( 14 ) 


2=1 


which is the well known AdS 2 x 


pd—1 


metric at zero temperature. 
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III. HOLOGRAPHIC FERMIONIC IMPLEMENTATION 


A. The flow equation 


In order to study the Fermionic theory dual to this bulk theory, we consider a probe spinor 
held with the charge q and mass m, which is dual to fermionic operator O in boundary CFTd 
with conformal dimension A = m + |. The action of spinor held is given by 

Sp = J (15) 

where F" = (e^)“F^ and the covariant derivative Da = da + — iqiAa with F^*^ = 

i[FM^F^], the spin connection {uj^y)a = a{&u)bi where (e^)“ forms a set of orthogonal 

normal vector bases. From the above Eq. ffTSj) . the equation of motion of Dirac held can be 
written as 

V^Da^! — rmp = 0 . (16) 


In order to investigate the Dirac equation in the Fourier space, we use the translational 
invariance i/j = J ijj. Then, we can rescale the Dirac held by "0 = 

and remove the spin connection completely in the Dirac equation. Considering the rotation 
symmetry in the spatial directions, we can simply set ki = kx- Then, it is convenient to 
consider the following basis 



where 0 / are two-component spinors, 1 is an identity matrix of size 2 ^* 2 ^ for odd d (or 2 2 
for even d). By taking the above operations, the Dirac equation ffT 6 l) can be written as 



+ y/^rna^ (g) 


+ qAt)ia'^ 1 T ka^ 0 [ ](^8) 

\(p2 J / 


Furthermore, by setting 0/ = 


Hi \ 

and introducing the ratios = —, the equation of 

X ^ V 

motion can be decoupled and reduced to 

+ qAt) + (-l)^/c] qAt) - (-l)^/c]^j, (19) 

V drr V Qtt V Qtt 

The in-falling boundary condition at the horizon for a; 7 ^ 0 is 


* r->-ro . 


( 20 ) 
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B. Near boundary and boundary conditions 


Near boundary r —)■ oo, the metric has the behavior as 

gtt Qrr , Qii r^. 

Substituting the near boundary metric into Eq. (flSD . we have 

\r‘^dr + rma^ — (cu + — (— (j)i = 0. 

Therefore, the Dirac equation becomes 

(r'^dr + rma^^ 0 / = 0 . 

The solution of the Dirac equation near the AdS boundary is 

where o/ and bj can be regarded as the response and the source, respectively, 
suppose the source and the response in Eq. fl2T|) are related by 



CLj 






the boundary Greens functions G{u, k) is given by 


M 


G = -*57°. 


Recalling the dehnition oi = ^, we can write the Green function as 


o I 6 0 

G{uj, k) = lim ' 

r^oo \ Q 


and the behavior of near the boundary is 


e* ™ r-^^Gii. 


( 21 ) 

( 22 ) 

(23) 

(24) 
Then, we 

(25) 

(26) 

(27) 

(28) 


7 





C. Spinor field at the near-horizon geometry and the dispersion relation 


Near the horizon region, the geometry approaches AdS 2 x ^ which is controlled by 
Eq. flT^ . Substituting Eq. flT^ into Eq. flTS]) . the Dirac equation Eq. flTHD can be expressed 


as 


dcMO = ^ 


Co 


c vco-c 


ma^(l)i{Q — i 


Co 


Co-C 


^q{d - 2)Ll f Co 




Co-C 


Co 


r* C VCo - C 

Near the AdS 2 boundary C —t 0, the above equation can be simplified into 

mL2 


CdMO = -uMO, u = 


rnjL2 - 


mjL2 + 


-TTiLq 


(^V/(C) 

aV/(C)- (29) 

(30) 


with TTij = —{—lyy. Following the discussion in jd], equation fl30|) is nothing but the 
equation of motion for spinor fields in the AdS 2 background with masses [m, rhj] with mj{I = 
1, 2) the time-reversal violating mass terms. In this case, the spinor field is dual to a spinor 
operators in the IR CFTi whose conformal dimensions of the operator are S = Uk + \ with 


= \ m'^Lo + mfEn — 


^l.q(d-2)Lly 

For the real Uk, the scaling exponent z of dispersion relation of our UV Green function 

6 | 

1 with Ukp > I 

^ with < i 


ImG 22 can be determined by 


X = 


(32) 


where Ukp can be obtained from Eq. fl^ at k = kp. The scaling exponent 2 ; contains 
important property of the Green function near the Fermi momentum k± = k — kp and it 
describes the excitation of dispersion relation 


uykp) ~ ky 


(33) 


where u^kp) is the real part of the complex frequency at which the Green function has 
a polej^. After determining the Fermi momentum from numerical calculation, we can 
analytically compute the scaling exponent z of the dispersion relation through Eq. fITT]) and 
Eq. fl32p . Note that z = 1 is an necessary character for the Fermi liquid. So in the following 
discussion, moves from 1 to other values should be the distinguishing feature for the phase 
transition from Fermi liquid to non-Fermi liquid. 




















FIG. 1: The 3D plot of Green function with a = O(left), a = l(middle) and a = 2(right), 
respectively. 


IV. FERMI SURFACE AND PHASE DIAGRAM 


In this section, we shall study the Fermi surface structure and the phase diagram (a,T). 
We will demonstrate that the transition from the Fermi liquid to the non-Fermi liquid 
happens when the impurity parameter a becomes larger. In numerical calculation, we only 
focus on the Green’s function G 22 (<^, k) due to the symmetry G 22 (<^, k) = Gii(a;, —fc), which 
can be easily analyzed from Eq. (fT9|) and the boundary condition Eq. fl20D . We set m = 0, 
g = 1 in most of our study unless we emphasize different setting. In addition, in the following 
discussion, we can set ro = 1 by rescaling. 


A. Results at zero temperature 

In this subsection, we will study the Fermi surface structure at zero temperature. We 
discuss the case with d = 3 as the hrst step. From Fig. [H we can see that a sharp quasi 
particle like peak emerges near w = 0 for different a, which denotes the Fermi surface. 
By dehnition, we can determine the Fermi momentum kp hy locating the position of the 
maximal value of ImG 22 near a; = 0 as shown in Fig. [2l We list some Fermi momentum kp 
with different a in Table [B We can see that with the increase of a, the Fermi momentum 
kp becomes smaller. After determining numerically the Fermi momentum kp, the scaling 
exponent .2 of the dispersion relation can be calculated by Eq. (j32ll analytically. The results 
have also been summarized in Table [T1 

The effect of the dimension d is explicitly reflected in the flow equation fiM]) . With careful 
numerical calculation, we list the results with different dimensions in Table [B For the same a, 
Fermi momentum is suppressed in higher dimensional dual theory. In addition, the scaling 
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ImG22 
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FIG. 2; Green function at tiny frequency with different a. 


d -1- 1 

4 

5 

6 

a 

0 

1 

2 

0 

1 

2 

0 

1 

2 

kp 

2.27694628 

2.06223335 

1.2532708 

1.8880638 

1.7800883 

1.3981258 

1.7228794 

1.6551679 

1.4260239 

Vk 

0.728524 

0.665326 

0.377719 

0.361108 

0.338395 

0.238421 

0.220036 

0.207121 

0.151011 

z 

1 

1 

1.32373 

1.38463 

1.47756 

2.09713 

2.27236 

2.41405 

3.31103 


TABLE I: kp-, Vk and z with different dimension and a (T = 0). 


exponent z of the dispersion relation decreases with the decrease of the spacetime dimension, 


which agrees well with the observation in Charged Gauss-Bonnet gravity in 


io|. 


From Table HI an interesting result is that the transition from Fermi liquid phase to the 
non-Fermi liquid phase happens when the impurity a is beyond some critical value. This 
motivates us to turn on the temperature and draw the phase transition diagram of a — T. 


B. Results at finite temperature 

Now, we turn to explore the Fermi surface structure at the hnite temperature and plot the 
phase diagram [a, T). Due to the thermal fluctuations at hnite temperature, the height of 
the quasi-particle peak becomes much lower and the width much broader (Fig. [3]). Using the 
same numerical method as the case at zero temperature, we present the Fermi momentum 
kp as the function of T with different a in Fig. 01 We can see that with the increase of the 
temperature, the Fermi momentum kp decreases for all a. 

Furthermore, we show the phase diagram {a,T) for d = 3 in Fig. [5l In Fig. [5l the line 
denotes the marginal Fermi liquid phase and the left region is the Fermi liquid while the right 
region corresponds to the non-Fermi liquid. It exhibits a phase transition from Fermi liquid 
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FIG. 3: Green function at different temperatures with a = 2. 



kf 


0.000 0.002 0.004 0.006 0.008 0.010 

T 


FIG. 4: Fermi momentum versus temperature for a = O(Blue), a = l(Red) and a = 2(Green), 
respectively. 


to non-Fermi liqnid indnced by the effective impurity a, which is qualitati vely consistent 
with the picture that a coherent/incoherent metal transition induced by a in 3l|. We would 
like to point out that the phase diagram {a,T) exhibiting in Fig. [5] seems to be opposite 
to the one relating to the doping drawn in the high-Tc cuprate superconductors ^ . This 
phenomenon somehow supports that the impurity a introduced here may not be considered 
as doping which is claimed in 43|] At least, a provides a new mechanism introducing 
impurity from holography. It would certainly be interesting if such a phase diagram can be 
observed in some real materials in the future. 

As is known that the model parameters m and q have effects on the types of dual liquids 


and t 
terns 


re phase transition from Fermi to non-Fermi liquids in the holographic fermionic sys- 


18| . It was addressed that with fermion charge q increasing, the Fermi momentum 


increases approximately linearly and the scaling exponent z of dispersion relation decreases 


43| treated a in a more general model as disorder-strength but the existence of 


^ The authors of 
Anderson localization needs to be further carefully confirmed. 
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FIG. 5: Phase diagram with d = 3. 


rapidly. Thus, the dual liquid changes from the Landau Fermi liquid type to the non-Fermi 
liquid is more difficult for larger charge. Here, we hud the similar tendency in our model. 
In the left plot of Fig. [6l we see that in the case with larger charge, the phase transition 
happens with stronger impurity parameter a. The reason is why larger charge is harder 
to make the phase transition happen so that q suppresses the effect of impurity. We move 
on to discuss the effect of the mass of fermions. For a given a, the Fermi momentum kp 
decreases and the scaling exponent increase with the increase of m which agrees well with 


the results in 


18| . At certain a, we always hnd a phase transition from the Fermi liquid to 


the non-Fermi liquid with hxed mass. From the right plot of Fig. [6l we see that at a hxed 
temperature, the critical a is smaller for larger mass, namely that larger mass promotes the 
phase transition. We also study the phase transition in higher dimension which may affect 
the type of the dual liquid. An interesting result is that the Fermi liquid phase is forbidden 
for > 0 in the d = 4 and d = 5 geometry. And it can be realized when we choose 
a negative enough which enhances the chemical potential of the system with a hxed 
temperature. However, this means that a should be imaginary, which should be further 
understood from holography. 


V. DISCUSSION AND CONCLUSIONS 

In this paper, we have explored the Fermi surface and the phase diagram (a,T) in the 
holographic framework based on the free probe fermions in Einstein-Maxwell gravity with 
spacial dependent massless scalar held, which introduces a new impurity ehects from holog¬ 
raphy in the dual held theory. The Fermi momentum was suppressed by the larger impurity 
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a a 


FIG. 6: Left: Phase diagram for d = 3 and m = 0 with different fermion charge. The line denotes 
Marginal Fermi Liquid(MFL) where z/fc = ^ for g = 0.9(Blue), q = l(Red), q = l.l(Green). Right: 
Phase diagram for d = 3 and q = 1 with different fermion mass . The line denotes MFL for 
m = O(Blue), m = 0.01 (Red), m = 0.02(Green). In each case, the left region of each line denotes 
Fermi liquid while the right side means non-Fermi liquid. 

parameter a while the conformal dimension of the Dirac held decreases as the impurity a 
increases at any temperature. Furthermore, we have demonstrated that the effective impu¬ 
rity parameter a can induce the phase transition from the Fermi liquid to the Non-Fermi 
liquid. This phenomenon results from that large positive a corresponds to lower chemical 
potential in the dual system with hxed temperature. 

Our work is a hrst step in studying the impurity effects on the fermionic spectral function 
in holography, in which the inhuence of the effective impurity or the spatial broken trans¬ 
lation symmetry only rehects on the Dirac equation through the metric correction. To gain 
more insights into the impurity effects of the fermionic spectral function in holography, we 
can further introduce an interaction between the massless scalar held and the Dirac held. 
This can produce rich Fermi surface structure, which can be compared with the results 
of Angle-resolved photoemission spectroscopy (ARPES) or Scanning tunneling microscopy 
(STM) in some real impurity materials, and more complete phase diagram including Mott 
insulating phase. The related work is under progress. It is also interesting to consider the 
non-relativistic Fermionic hxed point to see the effect of the impurity a on the flat band. 
We will present our results on this topic elsewhere in the future. 
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